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Abstract 

The phonon dispersion relation of NiTi in the simple cubic B2 structure 
is computed using first-principles density-functional perturbation theory with 
pseudopotentials and a plane- wave basis set. Lattice instabilities are observed 
to occur across nearly the entire Brillouin zone, excluding three interpenetrat- 
ing tubes of stability along the (001) directions and small spheres of stability 
centered at R. The strongest instability is that of the doubly degenerate ]VLy 
mode. The atomic displacements of one of the eigenvectors of this mode 
generate a good approximation to the observed 519' ground-state structure. 
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I. INTRODUCTION 



The realization of the technological promise of active materials in general, and shape 
memory alloys in particular, requires significant progress in the fundamental understanding 
of their nature and behavior. The relative chemical and structural simplicity of the shape 
memory alloy NiTi has made it a popular subject of both experimental and theoretical 
study@. 

The monoclinic f?19' low-temperature structure of NiTi is related to the cubic B2 high- 
temperature structure by approximately rigid shifts of alternate (110) planes along the (110) 
direction, resulting in a lowering of the symmetry and a corresponding change in the shape 
of the unit cell!. Previous first-principles studies have focused on accurate prediction of the 
ground-state structural parameters by relaxing forces and stresses for a particular choice of 
space group, complementing the experimental structural determination datai'S. 

The distortion that relates the two structures has the symmetry of a single M-point nor- 
mal mode of the B2 structure. This suggests a parallel between the martensitic transforma- 
tion from the B2 to the i?19' structure in NiTi and the paraelectric-ferroelectric "soft-mode" 
transition in perovskite oxides such as PbTiC>3, BaTiC>3 and KNbOsi. In these perovskites, 
first-principles calculation of the eigenfrequencies and eigenvectors of the unstable phonons of 
the cubic high-symmetry prototype structurei~i has shown that an excellent approximation 
to the ferroelectric ground-state structure is obtained by "freezing in" the single dominant 
unstable phonon (a polar optic mode at the zone center) with an accompanying change in 
the shape of the unit cell. The q-dependence of the instability across the Brillouin zone has 
been shown to determine the nature of local structural fluctuations in the high-temperature 
cubic paraelectric phasei and, through a first-principles effective Hamiltonian analysis, the 
transition temperature and related quantities!. 

In this paper, we take the first step in extending this approach to shape memory alloys 
by computing the phonon eigenfrequencies and eigenvectors of the cubic high-symmetry 
prototype structure of NiTi from first principles, as described in Section II. In Section 
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Ill, the dominant unstable phonon is identified and the structure it generates is compared 
with available knowledge of the ground state structure. We then investigate how these 
instabilities extend across the Brillouin zone. We conclude by speculating on the implications 
of this observation on the structure and properties of the high-temperature phase and on 
the martensitic transition. 



II. COMPUTATIONAL METHOD 

First-principles calculations of the structural energetics of B2 NiTi were carried out 
within density-functional theory with a plane-wave pseudopotential approach. Phonon 
eigenfrequencies and eigenvectors throughout the Brillouin zone were obtained using the 
Green's function formulation of density-functional perturbation theory (DFPT)i. The calcu- 
lations were performed with the PWSCF and PHONON co des§, using the Perdew-Zuneeio 
parametrization of the local-density approximation (LDA). Ultrasoft pseudopotentialslll 
for Ti and Ni were generated according to a modified Rappe-Rabe-Kaxiras-Joannopoulos 
(RRKJ) scheme^ with three Bessel functions^. The electronic wave functions were repre- 
sented in a plane-wave basis set with a kinetic energy cutoff of 30 Ry. The augmentation 
charges were expanded up to 660 Ry. The Brillouin zone (BZ) integrations were carried out 
by the Hermite-Gaussian smearing techniquell! using a 56 k-point mesh (corresponding to 
12x12x12 regular divisions along the k x , k y and k z axes) in the irreducible wedge. 
The value of the smearing parameter was cr=0.02 Ry. These parameters yield phonon fre- 
quencies converged within 5 cm -1 . The dynamical matrix was computed on a 6 x 6 x 6 
q-point mesh commensurate with the k-point mesh. The complete phonon dispersion rela- 
tion was obtained through the computation of real-space interatomic force constants within 
the corresponding boxS. 

In Table |, we report the equilibrium lattice parameter and elastic constants of B2 NiTi 
obtained from pseudopotential total energy calculations performed as described above, and 
from a previous pseudopotential calculation using a mixed basis set(MB)i. For comparison, 
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we also performed full-potential linearized-augmented-plane-wave calculations (FLAPW) 
within both the LDA and the generalized-gradient approximation (GGA)0. The lattice 
parameter and bulk moduli were obtained from a Birch function fit@. The discrepancy 
between the PWSCF and FLAPW values for the LDA lattice parameter is comparable 
(< 1%) to that between the two pseudopotential results. The FLAPW lattice parameter 
is slightly larger within GGA, as expected, and the bulk moduli show the typical decrease 
with increasing lattice parameter. Although these results are not directly comparable with 
the properties of the entropically-stabilized high-temperature B2 phase, characterized by 
large fluctuating atomic displacements, the experimental data included in Table [I] are in 
general agreement with the first-principles values. The DFPT calculations reported in the 
next section were performed at the PWSCF lattice parameter 5.594 a.u. 



III. RESULTS AND DISCUSSION 

The full phonon dispersion of B2 NiTi consists of six branches: three acoustic and three 
optic. Results of our calculations along the high-symmetry lines of the simple cubic BZ are 
shown in Figure |l[ The imaginary frequencies of the unstable modes are represented as 
negative values. The values of the phonon frequencies at the symmetry points are also listed 
in Table 0. All of the eigenvectors at the symmetry points T, X, M and R are uniquely 
determined by symmetry considerations except for the M5/ modes, to be discussed in more 
detail below. 

As expected on the basis of experimental observation of the low-temperature structure, 
the B2 structure is unstable at T = 0. The dominant instability is a doubly-degenerate 
M 5 / mode. Since there is a second M 5 / mode (at high frequency), the atomic displacements 
corresponding to the unstable M5/ mode are not uniquely determined by symmetry, but can 
be obtained only by diagonalizing the computed dynamical matrix. In addition, since the 
mode is doubly degenerate, there is a two-dimensional space of phonon eigenvectors that all 
correspond to the same frequency. Two particularly symmetrical eigenvectors in this space 
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are shown in Figure 0. In fact, the distorted structure produced by freezing the eigenvector 
in Figure 0(a) into the reference cubic structure yields, with an appropriate choice of overall 
amplitude, an excellent approximation to the observed ground-state £?19' structure. As can 
be seen in Table III, the main difference is in the angle 7, lowering the symmetry from 
orthorhombic to monoclinic. Energetically, this eigenvector is singled out from the others in 
this space not at quadratic order, but as the result of higher order terms, including strain 
coupling, not reported here. Detailed calculations will be presented in a future publication. 

While the evolution of the lattice instability away from the zone center cannot be directly 
obtained experimentally, Figure [l] shows that the unstable phonon branches actually extend 
throughout the Brillouin zone. From the M point, the instability extends more than half 
of the way to T and X. At T, the acoustic branches behave normally at very small q. As 
q increases along T-R and T-M, one and two individual optic modes, respectively, disperse 
strongly downward, mixing with the acoustic branches and then becoming unstable. As the 
R point is approached along T-R, the unstable longitudinal mode turns around and becomes 
stable once again. Along T-X the A± optic mode similarly disperses downward and mixes 
with the Ai acoustic branch, although along this line all phonon modes remain stable. 

This behavior can be better visualized in a three-dimensional view of the BZ showing 
the lo 2 =0 isosurface (Figure §). The region of stability is confined to three interpenetrating 
tubes along the cartesian axes, with additional bulges in the central region along the T-R 
directions, and to small approximately spherical regions around the R points. The unstable 
modes at M are thus continuously connected to the unstable mode along T-R. 

Phonon frequencies along symmetry lines in the high-temperature B2 phase have been 
experimentally measuredil and theoretically analyzed^. As explained in Section II, the 
properties of this phase are not directly comparable with calculations for the unstable zero- 
temperature B2 structure (note, however, that the calculated stable acoustic mode frequen- 
cies are in general agreement with their experimental counterparts). The high-temperature 
B2 structure is entropically stabilized and should have strong local distortions away from the 
high-symmetry average positions of the atoms. As in the previous work on KNbOs, the na- 
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ture of these disortions should be related to the distribution of unstable modes in the BZ, and 
could be experimentally characterized through diffuse X-ray scattering or with local probes, 
such as EXAFS. The soft-mode behavior of NiTi is also less straightforward than that of 
the ferroelectric perovskites. In the cubic phase, a soft TA mode at q — ^(§§0) freezes 
in to generate the intermediate-temperature R phase, and a softening of the M 5 / phonon 
associated with the ground-state B19' phase has not been directly observed. Theoretical 
investigation of the phase diagram and temperature-dependence of the phonon dispersion 
requires a first-principles effective Hamiltonian analysis, analogous to that for the perovskite 
oxides, beyond the scope of the present paper. 

IV. CONCLUSIONS 

In conclusion, we performed ab initio calculations of the phonon dispersion of the B2 
high-symmetry reference structure of NiTi. There are lattice instabilities throughout the 
entire Brillouin zone, with the dominant instability at M. A good approximation to the 
observed ground state structure of NiTi can be generated by freezing in a particular choice 
of eigenvector from the corresponding two-dimensional space. 
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FIGURES 

FIG. 1. Phonon dispersion for NiTi in the B2 structure with ao = 5.594 a.u. along symmetry 



lines in the simple cubic BZ. Symmetry labels are assigned according to the conventions of Ref. 21 
with Ni at the origin. The imaginary frequencies of the unstable modes are plotted as negative 
values. 

FIG. 2. Relative atomic displacements corresponding to two linearly independent eigenvectors 
of the unstable M' 5 mode at k = ^(330), transforming as (a) x + y and (b) y. Here, a portion 
of the B2 structure, with Ni atoms shown by filled circles and Ti atoms shown by open circles, is 
viewed along the (001) direction. All displacements lie in the x-y plane. 

FIG. 3. Zero- frequency isosurfaces of the phonon dispersion relation in the first octant of the 
Brillouin zone. Outside the interpenetrating tubes along the (001) directions and the small spheres 
centered at R. 
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TABLES 



TABLE I. 


Lattice parameter 


and elastic properties of NiTi in 


the B2 structure, ao in the 


parenthesis is 


GGA result. 










PW 


TV T T ~» / T ~"l P L \|\ 

MB(Ref. g) 


T — IT A r»TTTtl 

FLAPW a 


Exp. 


a (a.u.) 


5.594 


5.626 


5.561(5.670) 


5.698(Ref. [jJD 


B(Mbar) 


1.68 


1.56 


1.86 


1.40(Ref. |J) 


cn(dyn/cm 2 ) 


1.80 


1.68 




1.62(Ref. |J) 


Ci2(dyn/cm 2 ) 


1.68 


1.44 




1.29(Ref. H) 


C44(dyn/cm 2 ) 


0.77 


0.50 




0.35(Ref. 0) 



a This work, see Ref. 17 



TABLE II. Computed phonon frequencies of B2 NiTi with ao 
points in the simple cubic BZ. 



5.594 a.u. at the symmetry 



Label 


Frequency (cm 1 ) 


Label 


Frequency (cm 1 ) 


Lis 





Rl5 


30 


Lis 


200 


R25 


252 


X v 


77 


My 


91t 


X 5 / 


121 


Mr 


91 


x 5 


150 


M 4 / 


140 


Xi 


293 


Mr , 


264 
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TABLE III. Structural parameters of 52 NiTi distorted by freezing in of the eigenvector in 
Figure ||, a=2.960 A, 5=4.186 A, c=4.186 A, 7=90°, compared with the experimental monoclinic 
519' structure (Ref. ||), a=2.898 A, 6=4.646 A. c=4.108 A, 7=97.78°. While the distorted 52 
structure actually has space group Pmma, the structure is presented within P2i/m, the space 



group of the 519' structure, to facilitate comparison. 

Structure Wyckoff position x y z 

distorted 52 Ti(2e) 0.5 0.222 0.25 

Ni(2e) 0.678 0.25 

519' Ti(2e) 0.4176 0.2164 0.25 

Ni(2e) 0.0372 0.6752 0.25 
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